Polyelectrolyte-colloid complexes: polarizability and effective interaction 
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We theoretically study the polarizability and the interactions of neutral complexes consisting of a 
semi-flexible polyelectrolyte adsorbed onto an oppositely charged spherical colloid. In the systems 
we studied, the bending energy of the chain is small compared to the Coulomb energy and the chains 
are always adsorbed on the colloid. We observe that the polarizability is large for short chains and 
small electrical fields and shows a non-monotonic behavior with the chain length at fixed charge 
density. The polarizability has a maximum for a chain length equal to half of the circumference of 
the colloid. For long chains we recover the polarizability of a classical conducting sphere. For short 
chains, the existence of a permanent dipole moment of the complexes leads to a van der Waal's-type 
long-range attraction between them. This attractive interaction vanishes for long chains (i.e., larger 
than the colloidal size), where the permanent dipole moment is negligible. For short distances the 
complexes interact with a deep short-ranged attraction which is due to energetic bridging for short 
chains and entropic bridging for long chains. Exceeding a critical chain length eventually leads 
to a pure repulsion. This shows that the stabilization of colloidal suspensions by polyelectrolyte 
adsorption is strongly dependent on the chain size relative to the colloidal size: for long chains the 
suspensions are always stable (only repulsive forces between the particles), while for mid-sized and 
short chains there is attraction between the complexes and a salting-out can occur. 

PACS numbers: 82.70.-y, 82.35.Rs, 61.20.-p 



I. INTRODUCTION 

Polyelectrolyte(PE)-colloid complexes have recently 
motivated a great amount of computational and ana- 
lytical work jmjl & |_| J, @, § g 0, 0, H || 
U, H, HI, fL7f |§ ]l|ri§ H ■ The reason for this is 
not only the interesting behavior seen in such systems, 
but also their potential applications. For instance, some 
water-soluble paints are formed by PE-colloid complexes, 
in which a colloidal suspension is stabilized by adsorbed 
charged polymers. A better understanding of these com- 
plexes can give hints on how to avoid the precipitation 
of the particles from the suspension. Another example 
is the DNA-histone complexg, |^, ||, whose behavior is 
thought to be one of the crucial factors in the packing of 
DNA in living cells Q. 

Most of the previous work on PE-colloid complexes 
looked essentially at the structure and complexation 
behavior of the aggregates, e.g., of one PE and one 
macroion, g g 1 ff g |, || Q, Jj| OL M || M 0, 
one PE and many macroions |1C| , |12|, |15, , one 
macroion and many PEs [13|, or the effect of overcharg- 
ing of the macroion @, |9[ 16, [l8], |l9|. While some works 
discuss the interactions between two macroions in pres- 
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ence of either one or more than two short PE-chains 
[fill [20| , pl| we concentrate here on the complex formed 
by one spherical colloid and one PE as a whole and look 
at its polarizability, as well as the effective interaction 
p2[ between two isolated complexes. To do this we per- 
form Monte Carlo (MC) computer simulations of one and 
two neutral symmetric PE-colloid complexes, i.e., sys- 
tems where the charge of the PE exactly neutralizes the 
charge of the colloidal particle. We focus on the case 
where the bending energy is small compared to the en- 
ergy of Coulomb coupling, i.e., the chain is always com- 
pletely adsorbed on the macroion. 

While the results we will show are, strictly speaking, 
only valid in the limit of low density of complexes and 
zero salt concentration, we expect these to remain qual- 
itatively correct under more realistic conditions. As we 
observe in the simulations, the effect of salt at low con- 
centrations (put into the simulations through a Debye- 
Hiickel interaction) affects only slightly the results given 
for the salt-free system. Also, while both the polarizabil- 
ity and the interaction between complexes is very much 
influenced by the "microscopies" (i.e., by the structure 
of the chain), we will show that most of the computer 
simulation results can be understood using fairly simple 
principles. 

The paper is organized as follows: in section II we de- 
scribe the model used and the computer simulation de- 
tails. In section HI we present the simulation results for 
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PE-colloid complex and suggest analytical expressions for 
the polarizability for all chain lengths. Section IV deals 



with the interaction of two of these complexes for differ- 
ent distances and chain lengths. We describe how the 
force is calculated and the results are discussed. We pro- 
vide analytical results for the force for short chain lengths 
and long distances. Finally in section |y| we make some 
final remarks. 



II. THE SIMULATION MODEL 

We performed monomer-resolved Monte Carlo (MC) 
simulations using the model of Kremer et al. for single 
polyelectrolyte (PE) chains [f|| |IJ |||. The polyelec- 
trolyte chains are modeled as bead-spring chains of N 
Lennard- Jones (LJ) particles. This method was first ap- 
plied to neutral linear polymers and to a single star poly- 
mer [^4| later also for neutral |26|, [27j] and charged 
stars polymer systems [^8| ^9|. For good solvent con- 
ditions, a shifted LJ potential is used to describe the 
purely repulsive excluded volume interaction between all 
N monomers: 



^Lj(r) 




for r < 2 1 / 6 cr; 
for r > 2 1 / 6 cr. 



(1) 



Here, r is the distance of the interacting beads, a is the 
microscopic length scale of the beads and £lj sets the 
energy scale. We will from now on rescale all lengths 
with a according to 



~ r 
a 
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In accordance with previous work|26[ [28| , |29j| , we have 
chosen for the temperature T — 1.2elj/A:b, where k B is 
the Boltzmann constant. 

The connectivity of the bonded monomers is assured 
by a finite extension nonlinear clastic (FENE) potential 
acting between neighboring beads: 



V F (r) 



^) 2 ln 



1 



for r < r n 
for r > r„ 



(3) 



where &f denotes the spring constant and is set to /cp = 
7.0£lj- This interaction diverges at r = r m , which deter- 
mines the maximal relative displacement of two neigh- 
boring beads. The energy £lj is the same as in Eq. (Fy), 
whereas for the length scale r m we have chosen the value 
r m = 2.0. 

Finally, the full Coulomb interaction Vc (r) between all 
charged monomers has to be taken into account: 



Z-l 2j4 6^ 

V c (r ij ,Z i ,Z j ) = 3 



k B T£ B — , (4) 



where Zi and Z~ are the number of elementary charges 
of the interacting particles i and j, e is the elementary 



charge and ry the distance between the particles. In 
our polymer model, every monomer has charge valence 
z, meaning that the total charge valence of one PE is 
simply given by Z = N z. The Bjerrum length £b is 
defined as the length at which the electrostatic energy 
equals the thermal energy, viz. 



(5) 
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where eoe r is the permitivity of the solvent. For water in 
room temperature one obtains £b = 7.1 A. The solvent is 
only taken into account via the dielectric constant eoer- 
In our simulations, the Bjerrum length is fixed to £b = 
3.0, which is a realistic value for typical polyelectrolytes, 
such as the hydrophobic sodium poly(styrene-co-styrene 
sulfonate) (NaPSS) or the hydrophilic poly(acrylamide- 
co-sodium-2-acrylamido-2-methylpropane-sulfonate) . |3^] 
The interaction between one PE and a spherical 
macroion of radius i?o is modeled as follows: a monomer 
in distance r from the center of the colloid has a repulsive 
interaction V£j(r) of the truncated and shifted Lennard- 
Jones type 
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for r < Rq; 
for r > R , 



(0) 



and an attractive Coulomb interaction Vc(r, —Z, z), i.e., 
electrostatic interaction between a monomer of charge 
valence z and an opposite charge at the center of the 
macroion with valence Z — N z, which exactly neutralize 
the charge of the polyelectrolyte. Naturally, the repulsive 
interaction between two macroions each with charge va- 
lence Z in distance r is given by Vc(r, Z, Z). At a later 
stage we will include the effect of added salt by using a 
Debye-Hiickel (screened) interaction, via 
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ZiZj exp [—K(rij — Ri — Rj 

' /fry (1 + KRi)(l + KRj) 



(7) 



where Ri and Rj are the hard core radii of the two in- 
teracting particles, in particular Rq for the macroions 
and zero for the monomers and n is the inverse screen- 
ing length pi] ]. The inverse temperature is denoted by 
ft = 1/fceT. k is related to the salt molar concentration 



c s = K 2 /2N A h 



(8) 



where Na is Avogadro's number. We stress the fact that 
in the absence of salt, the system contains no small coun- 
terfoils: the chain and colloidal particle exactly cancel 
each others charge. 

In this work, we focus our attention on the case where 
the bending energy of the chain is small compared to 
the Coulomb energy of the complex. A good estimate of 
the bending energy of a chain with persistence length l p 
(isolated chain) and contour length I wrapped around a 
macroion to build a complex with radius R is given by 



E h 



end 



(9) 
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This should be small compared to the Coulomb energy, 
what yields 



2Z 



2*B 
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> 1 



(10) 



as the condition for a complete adsorption of the chain. 
We measured the persistence length of an isolated chain 
by calculating the bond vector correlations in dependence 
of the contour distance between the bonds and determin- 
ing the typical decay length. We find a linear dependence 
lp ps 1/2 for a wide range of contour lengths if no intrin- 
sic stiffness is applied. Inserting this scaling relation in 
( ^p| ) and using Z m Iz gives us a simple estimate for the 
wrapping state for chains with no intrinsic stiffness 



Az 2 l B R > 1. 



(11) 



Another limit where desorption can take place if the 
Coulomb coupling is small compared to the thermal fluc- 
tuations, e.g., the macroion is very big compared to the 
chain size. There is a complete adsorption of the chain 
onto the colloid if the electrostatic attraction is larger 
than the thermal fluctuations 



R 



> 1. 



(12) 



Strictly speaking, the adsorbed state is metastable: if the 
colloid and the finite chain are initially put far apart from 
each other, they interact essentially like two point-like 
charges. It is well known that the 1/r potential between 
two charges is not enough to bind them together, i.e., 
the Coulomb potential between two point charges is not 
enough to compensate for the loss of entropy through 
binding (cf. Manning condensation fl32]| ) . However, if the 
particles are put enough close to each other, the system 
needs a very long time before the unbound stable state 
is reached. To see this, let us assume that two oppositely 
charged particles are confined to a "spherical" box of 
radius L. The probability of finding them at distance r 
from each other is given by 



V(r) = 



exp{Z 2 £ B /r} 



(13) 



where Z = J R drV(r) and R is the contact distance be- 
tween the charges. The ratio between the probability of 
finding the particles at r' = r + £ B and at r is then given 
by 



eM-Z 2 l B /(r(r + e B ))} 



(14) 



For r = £ B and Z = 10, this corresponds to Rp ~ 2 x 
10~ 22 . In other words, if the two charges are initially put 
i B apart, the probability of finding them at 2l B apart is 
much smaller than to find them at l B . Note that R-p 
increases with r: if the particles are initially put very far 
apart from each other (say r/£ B of the order of Z), the 
system will flow to the stable state much faster, since the 



ratio Rp become of order unity. In our simulations, and 
in agreement with previous works, the PE has always 
stayed at the close vicinity of the colloidal particle, since 
the parameters are chosen such that Eq. ( |Io[) and Eq. (|T^ ) 
are always satisfied. 

Increasing chain stiffness has great influence on polariz- 
ability and interaction of these complexes. Some possible 
effects will be discussed briefly by introducing a simple 
harmonic coupling between the bonds of one PE 



N-l 



Vbond = ^ fcang(«i ~ 0!qY 



(15) 



where an is the angle formed between the connecting vec- 
tors of bonds i and i — 1 and the bonds i and i + X. We 
choose an equilibrium angle ao — ft and examine the ef- 
fects of increasing chain stiffness by enhancing the bond- 
ing constant fc ang . Unless if explicitly mentioned, /c ang is 
always chosen to be zero. 

We applied the standard MC Metropolis algorithm |3j| 
where one MC step consists of a trial move of every 
monomer. The macroions were kept fixed. The accep- 
tance rate for one MC sweep was approximately 50 per- 
cent. For faster equilibration and better statistics we in- 
troduced full chain moves where the whole PE is rotated 
around the closest macroion. Details describing these 
methods can be found in ]34| , |35|. Our systems were 
simulated without any cell in an infinitely wide space 
where the full accurate Coulomb interaction without an 
cutoff is taken into account. We observed that only at 
low Coulomb coupling or large intrinsic stiffness the ther- 
mal fluctuations lead to an entropic uncoupling of chain 
and colloid. After a long equilibration time (typically of 
1.5 x 10 5 to 2 x 10 5 MC steps), the different observables 
were averaged over runs with 10 6 to 10 7 MC steps. 



III. POL ARIZ ABILITY OF A SINGLE 
COMPLEX 

In this section we study the polarizability of a single 
complex consisting of one PE-chain with contour length I 
and linear charge density r = Ze/l [Z — N z and e is the 
elementary charge unit) which is adsorbed onto a spheri- 
cal colloid with radius Rq and charge valence Z. In order 
to calculate the contour length in the simulation, we mea- 
sure the mean distance between neighboring monomers: 
for a monomer charge z = 1 we obtain I ps Na and for 
z = 2 we obtain I « l.SNa, i.e., to = t = 1 for z = 1 and 
t = 1.54 for z = 2. In what follows we will extensively 
use as a parameter the ratio between the PE length and 
the effective diameter of the complex (2 R) 



I 

2R' 



(16) 



For x = 7r the chain length is equal to the circumfer- 
ence of the complex. For a chain much thinner than the 
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FIG. 1: Density profile of the monomers of one chain 
adsorbed to a spherical colloid with radius Rq = 5 
for different chain length x and different number of 
charges per monomer z plotted versus the distance 
r to the center of the colloid. With increasing x 
and z the monomers move closer to the center and 
the effective complex radius R decreases. See Fig. 
@ for simulation snapshots of one isolated complex. 




FIG. 2: Computer simulation snapshots of a PE adsorbed 
to a spherical colloid with radius Rq = 5 for (a) x = 0.79, 
(b) x = 2.78, (c) x = 8.20 and (d) x = 16.40. The 
number of elementary charges per monomer is z = 1. The 
monomers are rendered as dark spheres with diameter a. 



TABLE I: Parameters of the simulated system for z = 1 
and Rq = 5. We simulated chains with different monomer 
numbers N and chose a fixed complex radius R accord- 
ing to the density profiles, see Fig. x = 1/2R is 
the according chain lengths, while I = Na for z = 1. 
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colloid radius and at low temperatures, the effective ra- 
dius R should be given by R w Rq. However, since in 
our simulations the monomers have a size that is com- 
parable to the colloidal particle, this effective radius is 
approximately R « Rq + a, with some dependence on 
the number of monomers in a chain and on the charge of 
each monomer, as can be seen in Fig. [l]. In this figure, we 
show the density profile p(r) of the monomers surround- 
ing a colloid with radius Rq — 5 for chains with different 
number of monomers ./V and different monomer charges 
z = 1,2. As one would expect, a higher Coulomb cou- 
pling between the PE and the colloid (through a higher 
total charge in the PE) leads to a reduction of the ef- 
fective radius of the complex, here assumed to be at the 
maximum of the monomer distribution. Simulation snap- 
shots are shown in Fig. |^ for systems with different num- 
bers of monomers and no external field. As already men- 
tioned, the dependence of R on the parameters of the PE 
is so weak that it can be assumed for all practical pur- 
poses to be a constant for a wide range of chain lengths. 
In the following we will choose R — 6.3 for short chains 
x < tt and R = 6.1 for long chains x > tt as effective 
complex radius, for a fixed colloid radius Rq — 5.0 and 
2 = 1. See Tab. [|for a summary of the relation between 
complex radius, monomer number, and chain length for 
the simulated system with z = 1. The lengths for z — 2 
can be calculated accordingly; here we used R — 6.2 for 
x < tt and R = 6.0 for x > tt. The colloid valence is sim- 
ply given by Z — Nz. We checked that all forthcoming 
results obtained from the simulations are independent of 
above particularly chosen radius of the colloid. All the- 
oretical results in this and the following section should 
hold generally for all colloid sizes. 

We start our analysis with chains which are short com- 
pared to the size of the colloid, viz. x < tt. In this regime 
the PE-colloid complex has always a nonvanishing instan- 
taneous dipole moment, since the center of charge of the 
chain cannot coincide with the center of the colloid. We 
define the dipole moment P of the complex as 



N 



(17) 



where R c is the vector pointing to the center of the colloid 
and fj is the position of monomer i. Notice that, with 
this definition, the magnitude of P is given by Ze times 
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the distance between the center of the colloid and the 
center of charge of the PE-chain. It is easy to show that 
the mean dipole moment of a thermally fluctuating rigid 
dipole with length d and charge Ze in an electrical field 
E is 



P = Zed[ coth(/3 Zed E) 



1 



pZedE 



(18) 



where P is the magnitude of P and (3 = X/ksT, Intro- 
ducing the rescaled dimensionless dipole moment 



P* 



P 

ae 



and electric field 

E* = f3ae E 
one can rewrite Eq. (|l^) as 

P* = Zd(coth(ZdE*) - 



1 



ZdE* 



(19) 



(20) 



(21) 



where d = d/a. For low temperatures or strong electrical 
fields ZdE* ^ 1, the mean dipole moment saturates to 
P* = Zd, while for small fields ZdE* <C 1 and we recover 
the well known result 



2 J2 



p* _ 
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(22) 



For the case of a chain adsorbed onto a spherical colloid, 
the center of charge of the chain is on average somewhere 
between the center of the colloid and its surface. Assum- 
ing that the chain can be approximately described by a 
circular arc fluctuating on the surface of the colloid (see 
eg. Fig. ||(b)), the distance between the center of the col- 
loid and the center of charge of the chain is given by 



d = R 



sin(x) 



(23) 



where R = R/a. Our model is assumed to work for 
chains shorter than the circumference of the sphere, i.e 
Eq. (|3|) holds for < x < n. 

In Fig. [3] we compare the dipole moment of a complex 
in an external field for various chain lengths and z = 1 as 
a function of the external field from computer simulations 
and from Eq. (pi]). Note that Rq = 5 and the effective 
radius is assumed to be R = 6.3 for a; <■ it. For the 
theoretical effective dipole length d, we use the result 
from Eq. ( |2^ ) when < x < tt and Eq. [l(] is satisfied. If 
Eq. [h] is not fulfilled, the chain is essentially a straight 
rod that touches the sphere — in this case, the effective 
length of the dipole is simply given by R. On the other 
hand, for x > tt the chain cannot be approximated by a 
two-dimensional circle and Eq. ( pl| ) breaks down. 

The polarizability a of the complex is defined as the 
slope of its dipole moment as a function of the electric 
field 
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FIG. 3: Dimensionless dipole moment P* = P/cre of 
one PE-colloid complex versus the dimensionless exter- 
nal electrical field E* = f3ae E for different chain lengths, 
z = 1 and a colloidal radius Rq = 5. Symbols are 
computer simulation results, lines are theoretical results. 
Plot (a) shows the scaled chain lengths x = 0.79 (circles) 
and x — 1.19 (squares). The dipole moment increases 
with x. Lines are the theoretical results, Eq. (|2l|), us- 
ing d = R (solid) or the effective dipole length d — 
Rsm(x)/x (long dashed). Plot (b) shows the results for 
x = 1.59 (circles), x = 1.98 (squares), x = 2.38 (dia- 
monds), x — 2.78 (triangles). Note that the dipole mo- 
ment now decreases with increasing x. The long dashed 
lines is theory as used in (a). Plot (c) shows the results 
for x — 4.10(circles), x = 8.20 (squares), x = 16.40 (dia- 
monds), x — 24.60 (triangles). The solid line is the dipole 
moment of a conducting sphere with radius R cs = 6.1 ac- 
cording to P* — a cs E* , where a cs is given by Eq.(p6|). 
In (c) we also show the effect of increasing chain stiff- 
ness for x — 4.10. The harmonic coupling is quanti- 
fied by /3fc ang = 4.0 (crosses) and /3fc ang = 5.0 (plusses). 
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FIG. 5: Simulation snapshot of a PE-colloid complex 
in an external electrical field E. The chain has an in- 
trinsic stiffness due to the harmonic coupling Eq. (jig) . 
The bonding constant is /3fc ang = 4.0, the number of 
monomers is N — 100, E* = 0.4, z = 1, and R = 5.0. 



FIG. 4: Dimensionless polarizability a versus the 
scaled chain length x at small applied fields. The 
symbols denote the simulation, the dashed lines are 
the theoretical results from Eq. ( p5| ) (inset) and 
Eq. ( |26| ) for charges per monomer z — 1 (circles) and 
z = 2 (squares). The colloid radius is R = 5. 



The behavior of a is shown in Fig. [|: in the regime of 
small electric fields, the polarizability is constant, here 
we are in the linear response regime, see Eq. (^2). 

In the regime where the interaction energy between 
the dipole and the applied field largely exceeds the ther- 
mal energy (ZdE* 3> 1) the polarizability becomes very 
small; this latter behavior is not surprising, since in this 
regime the dipole is essentially pointing towards the field, 
and thermal fluctuations which tend to destroy this align- 
ment at small fields become unimportant. In our simula- 
tions we never really reach the point where the polariz- 
ability becomes zero since at some point the field is strong 
enough to induce the unbinding of the chain from the 
colloidal particle. Plot ||(a) shows the dipole moment for 
two short chains, viz. x — 0.79 and x = 1.19. Notice that 
the system with former value of x is very well described 
by P with d = R, while the latter tends to approach the 
dipole moment which uses d given by Eq. (p3|). The rea- 
son for this is the fact that in the system with smaller x 
the polyelectrolyte has a bending energy of the order of 
the Coulomb energy of the complex, see Eq. (|Io|), and 
the system behaves more like a rigid rod in contact with 
a sphere than as a circular arc, as described by Eq. (23). 
The curve for x = 0.79 in (a) stops abruptly since the 
critical field for the desorption is reached. 

Longer chains lead to a decrease of the polarizability at 
small fields as can be seen in plot (b) of Fig. 0, where we 
plot x = 1.59, 1.98, 2.38 and 2.78. This effect is a result 
of the dislocation of the center of charge of the chain to a 
position that is closer to the center of the colloid, accord- 



ing to Eq. (|23|). Note that complexes with medium-sized 
chains tend to approximately exhibit dipole moment of a 
circular arc on a sphere. However, complexes with long 
chains, x > ir tend to have a linear response to a wider 
range of the applied field (see plots (b) and (c)), which 
is a consequence of the fact that larger chains have more 
charge at fixed charge density, and so the fields needed 
to deviate from a linear law are larger; also, the fact that 
longer chains have more possible conformations leads to 
a qualitative change in the response of the complex, and 
a simple model as the one used above is no longer valid. 

Let us have a closer look at the polarizability for small 
fields and short chains. Inserting Eq. ( p3[ ) into (p2|), one 
obtains, for small fields, the dimensionless polarizability 
of the complex as a function of the PE-chain length and 
charge density 



4r 2 E 4 . 2 
a = sin (x). 



(25) 



where the linear charge density r has been rescaled such 
that f = err. For a fixed t, this polarizability shows 
a maximum at x = 7r/2 for a chain with length equal 
to half the circumference of the colloid and vanishes for 
x = ir, when a chain has a full circumference. Indeed, 
this behavior predicted by theory is very close to in the 
computer simulation results in the inset of Fig. |[ where 
we plot a for small fields for a colloid with fixed radius 
Ro = 5 and R = 6.3 and monomer charges z = 1 and 
z = 2. 

In the cases where the complex contains a long chain 



the polarizability Eq. (E3) loses its validity. Note 
that in this case, the chain is wrapped around the sphere 
more than one time and one expects a small or van- 
ishing dipole moment. Because the monomers tend to 
distribute themselves on average more or less homoge- 
neously throughout the surface of the sphere, one could 
expect the response of the PE-colloid complex to electric 
fields to become similar to the response of a neutral con- 
ducting sphere, which has a dimensionless polarizability 
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given by (see. e.g., p6| ) 

a cs = ^ (26) 
£b 

where R cs is the radius of the sphere and £b — ^b/c. Not 
surprisingly, in the limit of very large chains i > 7r this 
latter result is recovered, as shown in Fig. [|. For this 
comparison we have chosen R cs = 6.1 for the effective 
sphere radius. Note that the simulation results for long 
chains are independent of the linear charge density in 
agreement with Eq. (p6|). 

Let us now discuss possible effects of an increasing in- 
trinsic stiffness of the chain. For short chains we observe 
an increasing polarizability if the harmonic coupling is 
strong enough to decrease the curvature of the chain 
around the spherical macroion meaning that the bend- 
ing energy reaches values of the order of the Coulomb 
energy. . 

The center of charge of the PE is less close to the colloid 
center and the effective dipole length increases. As an 
upper limit the maximal dipole moment is then given by 
Eq. (^|) inserting d — R, with R effective complex radius. 
However, for the case, where the bending energy comes 
close to the Coulomb energy, we observe that the chain 
decouple easily from the colloid and no bound complex 
is formed. 

More interesting is the behavior for long chains. For 
long chains with a large bending energy and no electri- 
cal field, so called one-tail or two-tail configurations sets 
in§, |, y, [Si ||. Here part of the PE is still adsorbed 
on the colloid while one or two rigid parts are stretched 
away from the macroion to lower the bending energy and 
so minimize the total energy. In an electrical field the 
one tail configuration is obviously promoted due to the 
strong dipolar interaction with the external field, see Fig. 
U for a simulation snapshot. We make the following ob- 
servation in our simulation: if the chain is fully adsorbed 
on the colloid the polarizability is still given by Eq. ( |2^ ) 
independent of the higher chain stiffness, only the effec- 
tive complex radius increases slightly due to a less tight 
wrapping around the colloid. But for a sufficiently strong 
chain rigidity a small external electrical field can now pull 
a part of the chain in direction of the chain to lower the 
bending and dipole energy and a dipolar one tail config- 
uration is formed. Two examples of the dipole moment 
are shown in Fig. [3] (c) for x = 8.20, /3fc ang = 4.0 and 
(3kan B = 5.0. For small fields the chain is fully adsorbed 
and the polarizability is the same as for the /?fc ang = 
case. Then the increased stiffness leads to a continuous 
dewrapping in the one-tail configuration and the polariz- 
ability increases dramatically due to formed dipole. Fur- 
ther increasing of the electrical field lengthens the tail 
and eventually decoils the chain from the colloid. The 
longer the tail, the larger is the dipole moment and the 
polarizability increases with growing field until complete 
desorption is accomplished. A snapshot for an one-tail 
complex in an electrical field E is shown in Fig. ra. 
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FIG. 6: Simulation results of the dimensionless force 
F* = (3 a F between two PE-colloid complexes ver- 
sus the surface-surface distance ? for different chain 
lengths x and z = 1. The colloid radius is R = 
5.0. Lines are guide to the eye. Plot (a): short 
chains x % it. Plot (b): long chains x > ir. 



Finally, adding salt obviously increases the polarizabil- 
ity as the range and strength of the Coulomb attraction 
between PE and the colloid is screened and the effective 
complex radius grows. We take the presence of salt into 
account by using a Debye-Hiickel potential between the 
monomers (and ignoring any polarizability due to the 
salt), going up to a salt concentration of c s » 1.0M. We 
observe that Eq. (|2^) for short chains and Eq. (|2^) for 
long chains still hold as long a stable complex is formed 
but with a slightly larger effective radius R or R cs , re- 
spectively. 

IV. EFFECTIVE INTERACTION BETWEEN 
TWO COMPLEXES 

We now turn to the interaction between two complexes 
at distance r. The effective force acting on the center of 
one colloid in the simulation is calculated as follows: con- 
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(c) 



• # 

# 



by performing the statistical average 




FIG. 7: Snapshots of two interacting PE-colloid com- 
plexes with colloid radius Rq = 5, z = 1, chain length a; 
and surface-to-surface distance s. From top to bottom: 
(a) x = 0.79, 8 = 8.0; (b) x = 0.79, s = 1.5; (c) x = 1.59, 
s = 10.0; (d) x = 1.59, s= 1.5; (e) x = 4.10, s = 2.5; (f) 
x = 8.20, s = 2.0; (g) x = 3.27, s = 11.25, /% ang = 4.0 



sider two spherical colloids with radius Rq and center-to- 
center distance r, which means that they have a surface- 
to-surface distance of s = r — 2Rq. The energy of one col- 
loid in presence of one monomer at distance r m is the sum 
of Coulomb attraction and short-ranged LJ-repulsion 



M N 



F m (r) 




(28) 



where the sum in i goes over all M PE-chains and the 
sum in j considers all N monomers of each chain. r[™ is 
the distance of monomer j from the ith PE-chain to the 
center of the considered colloid. The brackets in Eq. j28| ) 
denote ensemble statistical average. We obtain the total 
force between the colloids in distance r by adding the 
direct repulsive Coulomb force between the colloidal par- 
ticles, viz. Z 2 e 2 /47reoe T .r 2 . The total force can be then 
divided into two contributions, 



F(r) = F c (r) + F LJ (r). 



where 



F c (r) = 



Z 2 e 2 
After 2 



M N 



\i=l j = l 



r ij: Z ) 



(29) 



(30) 



is the Coulomb force due to the electrostatic attraction 
of all monomers and the electrostatic repulsion of the 
second colloid and 



M N 



(31) 



= 1 3=1 



resulting from the excluded volume interaction of the 
monomers through the aforementioned Lennard- Jones 
potential. Since this force is essentially a hard-core ex- 
clusion force, one can consider this as describing the in- 
teraction due to gain and loss of entropy. For this reason 
we call it entropic force. 

Based on the previous discussion on the polarizabil- 
ity of a single complex, one can essentially expect two 
distinct regimes for the interaction between two equal 
symmetric complexes, viz. when the chains are short and 
when the chains are long. Fig. ^| summarizes our results, 
where the dimensionless force F* = (3aF is plotted for 
systematical increase of x and z = 1 as a function of s, 
the surface-to-surface distance. Negative forces mean at- 
traction while positive forces mean repulsion. See Fig. |^ 
for some typical simulation snapshots. Notice that for 
short chains, a long-range attractive tail is present. This 
is consistent, as we will demonstrate, with a dipolar in- 
teraction between the complexes. For longer chains, this 
long range tail disappears, and only a strong short-range 
attraction (also present for shorter chains) survives. In 
the limit of very long chains, only a repulsive steric inter- 
action is present. We discuss these results in more detail 
in the following. 



W m ) = Vc(r m , -Z, z) + V%{r m ). (27) 



Short chains 



After a sufficiently long equilibration time, the force F m 
of all monomers acting on one of the colloids is measured 



In the previous section we showed that the complexes 
with short chains respond like permanent dipoles to elec- 
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trie fields. When interacting with each other, it is rea- 
sonable to expect that in such case a dipole-dipole in- 
teraction shows up, at least at large separations between 
complexes. One way of estimating this long-range tail 
is by looking at the interaction between two fluctuating 
permanent dipoles P = Zed with length d as given by 
Eq. (E3) expressed through the so called Keesom energy 



„ , , 1 ZH\ d 4 



(32) 



After some algebra one obtains for the interaction energy 
between two PE-colloid complexes with < x < 7r 
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f3v(r) = — - 



■ sin 4 (a;). 



(33) 



where r is the distance between the dipole and the po- 
larizable object. Accordingly, the attractive force acting 
on the polarizable object is 



F*(r) 



32 U 



r 4 R 8 



sin 4 (a;). 



(34) 



FIG. 8: Interaction force between two PE-colloid 
complexes for long distances and short chains (x < 
7r). The plot shows a comparison made between 
the analytical expression Eq. ([mJ) for the dipole- 
dipole interaction (lines) and the computer simula- 
tion results (symbols) for colloids with size Rq = 5 
and various chains lengths x and monomer charge z. 



We easily see that the electrostatic interaction has a max- 
imum of strength for a length x = 7r/2 according to the 
behavior of the polarizability of one isolated maximum. 
Another remarkable feature is the strong dependence on 
the complex radius which goes with the eighth power for 
a fixed x. A comparison of expression ([m]) with the long 
distance tail of the force between the complexes is done 
in Fig. H for different chain lengths and charge densities, 
leading as expected to a nice agreement. As theoreti- 
cally predicted the simulations show a maximal dipole 
interaction for a length x ~ n/2. 

In order to obtain some insight into the conformational 
properties of the interacting complexes we introduce the 
order parameter 



P1P2 
\P1WP2 



(35) 



which is the normalized dot product between the two 
dipole moments of each complex. For a dipole-dipole in- 
teraction — 1 < <fr < 1 should be positive as both dipoles 
show in the same direction. Negative values of <fi mean 
that the dipoles are on average pointing in opposite di- 
rections. By studying the average distance between the 
center of charge of the chains in our simulations, we found 
that every time <fi becomes negative the monomers accu- 
mulate between the colloids, meaning that each dipole 
moment points in direction of the other complex in aver- 
age. 

Fig. H shows as a function of separation for different 
chain lengths with z = 1 at fixed complex radius. As ex- 
pected, for short chains and large distances, <f> is positive: 
the dipoles of the two complexes basically point in the 
same direction and produce the dipole-dipole interaction 




ililillmiiii 
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FIG. 9: Simulation results (symbols) of the order pa- 
rameter 4> versus the surface-to-surface distance of two 
PE-colloid complexes with Rq = 5.0, z = 1, and differ- 
ent chain lengths. The solid line is the numerical solu- 



tion of the toy model sketched in Fig. 10 for x = 0.79 



as discussed above. Note that the largest values for <f> 
are obtained for complexes with chains with x = 1.59, a 
value close to x — n/2 where the polarizability is max- 
imal. Typical configurations in the dipole-dipole regime 
are shown in Fig. ^ (a) and (c) where the order parameter 
4> has positive values. A relatively sharp transition from 
positive to negative values of <f> is observed for distances 
s w 4 for all values of x. 

For small distances, 4> becomes negative for all sys- 
tems studied. As mentioned before, this means that the 
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centers of charge of the chains are located somewhere be- 
tween the colloids, something that turns out to be always 
true as s < 4. This can also be seen in snapshots for close 
distances of the complexes as for instance in Figs. 0(b) 
and (d). 

In order to obtain a better understanding of the small 
distance behavior of the force for complexes with short 
chains, Fig. || (a), consider the toy model depicted in 
Fig. [h]. Here the chains are modeled as hard spheres with 
effective diameter a. Each sphere sticks to one colloid and 
can only move on the surface of the macro sphere with 
radius Ro- One small sphere belonging to one colloid can 
not penetrate into the other complex, and the center-to- 
center distance between the colloids is r. In this case, the 
effective complex radius is R = Ro + a/2. 

As previously shown, complexes with short chains be- 
have like rigid dipoles with length d which is dependent 
of the chain contour length I. Hence, in our toy model 
the complex has a dipole moment in direction from the 
center of the colloid to the center of the small sphere 
with an effective length d and charge Z. For distances 
r > 2(Ro + a) this model shows a simple dipole-dipole at- 
traction, whereas for smaller distances the dipole-dipole 
attraction is modified by the restricted number of avail- 
able configurations. 

One can find the force between the complexes for this 
toy model by numerically solving the partition sum and 
calculating the free energy as a function of the distance 
r. The results for a chain length x = 0.79 is plotted in 
Fig. [ll] together with the corresponding simulation data. 
According to the results from the previous section we use 
d = 6.3. We use the diameter of the small spheres a as 
a fit parameter, which should in any case be of the order 
of the monomer size. This parameter does not affect the 
dipole-dipole interaction, being only important when the 
deviations from the pure dipole-dipole interaction sets in. 
For x = 0.79 and z = 1 we used a = 3.6. This model 
leads to a qualitative agreement for the short range be- 
havior of the force, as depicted in Fig. [ll]: while the scale 
of the force is larger for the toy model, the minima and 
maxima are located at approximately the right distances. 
In other words, this means that the behavior observed for 
the force at small distances and for small chains (a force 
with two minima) is at least partially explained by the 
restriction of the number of available configurations of a 
system that is still essentially composed by two interact- 
ing dipoles. 

It is also interesting to study the behavior of the or- 
der parameter <f> of our toy model solved for the length 
x = 0.79 and z = 1 also plotted in Fig. ^. The qualita- 
tive behavior in dependence of the distance is quite the 
same as in the simulation. At s » 4.5 the dipole-dipole 
regime breaks down and <f> descends to negative values for 
closer distances s. The negative value of <f> means that 
in our toy model both small spheres are located mainly 
between the macroions. One concludes that due to the 
excluded volume interaction of the small sphere/chain 
with the macrosphere the dipole-dipole interaction is dis- 




FIG. 10: Schematic representation of the toy model 
for two PE-colloid complexes at distance r for short 
chain lengths (x < ir). One complex behaves 
like a dipole with charge Z and length d (dumb- 
bells). The short chain is modeled as a simple 
sphere sticking on the surface of the colloid in dis- 
tance R = Rq + a/2 where Ro is the colloid radius 
and a is the effective diameter of the fluctuating chain. 
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FIG. 11: Interaction force between two PE-colloid 
complexes with a short chain x = 0.79 and z = 
1. Simulation results (circles) are compared to 
the numerical solution of the toy model depicted 
in Fig. 10 (dashed line) and the analytical re- 
sult for the dipole-dipole interaction (dotted line). 



turbed and the the spheres/chains locate themselves be- 
tween the macroions, contrary to what happens in the 
dipole-dipole interaction regime. 



B. Long chains 

Increasing the chain length we observe a deep short 
ranged attraction as plotted in Fig. (b) which is in 
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FIG. 12: Plot (a): simulation results of the entropic part 
of the force F£j = /3aF LJ defined in Eq. @ for dif- 
ferent chain lengths and z = 1. The colloid radius is 
i?o = 5.0. Plot (b): simulation results for the Coulombic 
part of the force F£ = a F c Eq. ( |30| ) using the same 
parameters as in plot (a). Lines are guide to the eye. 



FIG. 13: Simulation results of the force between 
two colloids with R Q = 5.0 carrying two chains of 
length x = 4.10 (squares) and one chain with dou- 
bled length x — 8.20 (circles). The monomer va- 
lence is z = 1. Lines are guide to the eye. 
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agreement with previous studies) 20, |2 



|. The depth 

of attraction initially grows with increasing chain length 
while the range of attraction shrinks. Further increasing 
of the chain length leads to a pure repulsion of the com- 
plexes. The attraction is essentially of entropic origin, 
as the comparison between the Coulomb force Fc and 
the entropic force Flj as defined in Eq. ( |30| ) and j3l| ) 
shown in Fig. p"2| demonstrates: for x = 3.27 (large chain) 
the entropic force clearly dominates over the Coulomb. 
For distances comparable to some monomer lengths the 
chain is able to gain entropy by bridging the two colloids 
with some parts of the chain, as seen for instance in snap- 
shots Figs. |^(e) and (f). In agreement with the absence 
of the dipole-dipole interaction between the complexes, 
note that for chains with x > 7r, only zero or negative 
values of <f> are observed: at short enough distances, the 
chains bridge between the colloids, otherwise the com- 
plexes basically ignore each other and no long-range force 
is present. 



FIG. 14: Simulation results of the force between two 
PE-colloid complexes for a fixed chain length x — 
3.27, z = 1, and different rigidity k = f3k ang ac- 
cording to Eq. jl5|). Lines are guide to the eye. 



Not surprisingly, for small distances the attraction of 
two long chains is indistinguishable from the behavior of 
one chain with the combined length of the two chains, 
as depicted in Fig. [l3], where the force between two com- 
plexes with chains with x — 4.10 is compared with the 
force between to colloids in the presence of only one chain 
with x = 8.20. The force is the same for small distances 
s < 5 where the entropic bridging regime is found. The 
linear long ranged attraction of the single chain system 
is due to an energetic bridging already observed and dis- 
cussed in Refs. @ and 

For complexes with very long chains only a pure re- 
pulsion is observable, see eg. x = 32.79 in Fig. || (b). If 
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the chain is long enough to cover the whole surface of the 
macroion no bridging is possible due to the electrostatic 
repulsion between the approaching chain and the shell 
of monomers around the macroion. A simple theoretical 
estimate for the critical length when pure repulsion oc- 
curs is obtained by equating the available surface of the 
spherical complex A = 4wR 2 with the area A — la of 
the covering chain with length I and thickness a, which 
yields 

x c = 2nR/a. (36) 

We compared above estimate for two different complex 
radii R s=s 3.0 and R w 6.0 with simulation results and 
found x c ks 15 and x c ~ 33 which is very close to the 
theoretical results x c = 18.85 and x c — 37.70. See Fig. 
^ (b) for examples for the force between two complexes 
with a colloid radius Rq = 5. 



C. Effects of chain stiffness or salt 

Increasing chain stiffness has an enormous influence on 
the interaction for all chain lengths. The mechanism of 
interaction changes completely as now the rigid chains try 
to bridge between the colloids to lower their bending en- 
ergy which now is not negligible as compared to the elec- 
trostatic energy. Simulation results for increasing chains 
stiffness are shown in Fig. [l4| for a chain length x = 3.27. 
At a critical distance the chains bridge from one colloid 
to the other as shown in the snapshot in Fig. |?] (g). We 
checked that the isolated complex with this length and 
bonding parameters shows always a complete adsorption 
of the PE to the charged sphere and no one- or two-tail 
configurations. The functional behavior of the force is 
linear for a wide range of the distance. This is similar to 
the interaction of two colloids with only one chain, see 
Fig. |l^ discussed in Ref. [^0| where the sharing of the 
chain was named energetic bridging as one chain is then 
able to neutralize both colloids and lower the electrostatic 
energy. In our similar case with two chains the bridging is 
due to lowering of the bending energy of the rigid chains. 
Adding salt shows no surprising behavior of the inter- 
actions for all chain length. The shape and qualitative 
behavior of the force remains the same for low salt con- 
centration, only the magnitude of attraction decreases. 
As expected, the interaction of short chains is more af- 
fected by adding salt since the long-range attraction is 
due to electrostatic interaction which are now screened. 
Here differences already occur for a salt concentration of 
c s k, 0.05M and the long-ranged attraction vanishes com- 
pletely for c s w 0.2M (measured for x = 1.59 and z = 1). 
The interaction of long chains has entropic origin and re- 
mains unaffected till c s exceeds values of approximately 
0.1M for a chain length x = 4.10 and z = 1. For values 
of c s w 1.0M still a weak attraction is observed. 



V. FINAL REMARKS 

In summary, we studied the polarizability and interac- 
tion of PE-colloid complexes for the case where the chain 
is fully adsorbed on the macroion. The polarizability is 
strong for small chains and small electrical fields and has 
a maximum for chain length comparable to the circum- 
ference of the colloid. It can theoretically be described 
in terms of a thermally fluctuating dipoles with an ef- 
fective length which depends on the chain length. For 
chains longer than the circumference of the macroion the 
polarizability is of the order of the polarizability of a 
classical conducting sphere with radius of complex size. 
The interaction shows a van der Waals like attraction 
for short chains and a deep short ranged attraction for 
midsized chains. Exceeding a critical chain length x c 
the interaction is completely repulsive. One could in 
principle try to explain the strong short-range attraction 
seen in the simulations through the attractive interac- 
tion seen between equally charged walls in the presence 
of polyelectrolytes or counterionsfliol (using the Der- 
jaguin approximation to correct for the spherical geom- 
etry). However, one should note that the Derjaguin ap- 
proximation is only valid in the limit where the colloidal 
particle is much larger than the adsorbed chains, which 
is not the case in the simulations presented here. 

One of the interesting conclusions that one can draw 
from this work concerns the stability of colloidal sus- 
pensions which are stabilized by adsorption of polyelec- 
trolytes. While for chains which are very long compared 
to the size of the particles x > x c this is indeed an ef- 
fective mechanism of stabilization (as shown in Fig. ^ 
for x = 32.79, where the forces are always repulsive or 
zero), complexes with chains with a length comparable 
with the colloidal size exhibit a short range attraction 
due to the so-called bridging between the particles. How- 
ever, for chains that are smaller than the colloid, a long- 
range attraction is present and one should expect such 
suspensions to behave like dipolar fluids pLT|]. This shows 
that the stabilization of colloidal suspensions by polyelec- 
trolyte adsorption is strongly dependent on the chain size 
relative to the colloidal size: for long chains the suspen- 
sions are always stable (only repulsive forces between the 
particles), while for mid-sized and short chains there is 
attraction between the complexes and a salting-out can 
occur. Although our studies were essentially done for iso- 
lated, strongly charged and symmetric complexes these 
conclusions should hold in more realistic situations where 
the chains have some polydispersity and the complexes 
are at relatively high concentrations, since the only key 
point is the existence of the complexes. 
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